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We propose a phenomenological spin fluctuation theory for antiferromagnetic quantum tri- 
critical point (QTCP), where the first-order phase transition changes into the continuous one 
at zero temperature. Under magnetic fields, ferromagnetic quantum critical fluctuations de- 
velop around the antiferromagnetic QTCP in addition to antiferromagnetic ones, which is in 
sharp contrast with the conventional antiferromagnetic quantum critical point. For itinerant 
electron systems, we show that the temperature dependence of critical magnetic fluctuations 
around the QTCP are given as XQ ^ T^'^^'^ (xo oc T~^^'^) at the antiferromagnetic ordering 
(ferromagnetic) wave number q — Q (g = 0). The convex temperature dependence of Xo^ is 
the characteristic feature of the QTCP, which is never seen in the conventional spin fluctuation 
theory. We propose that the general theory of quantum tricriticality that has nothing to do 
with the speciflc Kondo physics itself, solves puzzles of quantum criticalities widely observed in 
heavy- fermion systems such as YbRh2Si2, CeRu2Si2, and /3-YbAlB4. For YbRh2Si2, our theory 
successfully reproduces quantitative behaviors of the experimental ferromagnetic susceptibil- 
ity and the magnetization curve by choosing the phenomenological parameters properly. The 
quantum tricriticality is also consistent with singularities of other physical properties such as 
speciflc heat, nuclear magnetic relaxation time 1/T\T, and Hall coefficient. For CeRu2Si2 and 
/3-YbAlB4, we point out that the quantum tricriticality is a possible origin of the anomalous 
diverging enhancement of the uniform susceptibility observed in these materials. 

KEYWORDS: quantum critical phenomena, self-consistent renormalization theory, tricritical point, quan- 
tum tricritical point, heavy-fermion systems, YbRh2Si2, CeRu2Si2, /9-YbAIB4 
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1. Introduction 

In strongly correlated electron systems, energy scales 
of interactions among electrons become comparable to 
those of band widths. Due to such strong correlations, it 
is often observed that several phases (for instance, nor- 
mal metals, magnetically ordered phase, and supercon- 
ducting phase) compete each other. As a consequence of 
this competition combined with quantum fluctuations, 
critical temperatures of phase transitions often become 
zero and a quantum critical point (QCP) emerges. ^""^ 

Heavy-fermion materials are suitable systems for the 
study of the QCP. In heavy-fermion systems, because of 
the competition between two energy scales,'* namely, the 
Kondo temperature and the Ruderman-Kittel-Kasuya- 
Yosida interaction, critical temperatures of magnetic or- 
der often become zero and the QCP appears. Actually, 
QCP has widely been observed in heavy-fermion materi- 
als by controlling pressures, external magnetic fields, and 
chemical substitutions.^""^ 

Near the QCP, it has been proposed that quantum 
fluctuations modify the electronic properties drastically. 
More concretely, electrons do not follow the fundamen- 
tal and textbook properties of Landau's Fermi liquid, 
which are universally seen in normal metals. This un- 
conventional behavior is called non-Fermi-liquid behav- 
ior The non- Fermi- liquid properties of metals have 
attracted much interest because novel quantum phases 



including exotic superconductors are found in the region 
where such non-Fermi-liquid behaviors are observed. 

Around thirty years ago, one of the standard picture 
to understand the non-Fermi-liquid properties was es- 
tablished,^ where non-Fermi-liquid properties have been 
successfully explained in various cases by spin fluctua- 
tions around the QCP of the ordinary second-order tran- 
sition in the framework of Ginzburg, Landau and Wil- 
son.*"^ We call this standard theory conventional spin 
fluctuation theory which covers so-called self-consistent 
renormalization (SCR) theory ^' ^ and rcnormalization- 
group treatment.^'* In Sec. 2, we will review how the 
conventional SCR theory describes the non-Fermi liquid 
behaviors near the QCP. 

However, it has been pointed out that this standard 
picture does not explain recent many experimental re- 
sults. Even when the apparent QCP is seen, physical 
properties do not follow the prediction of scalings by the 
conventional spin fluctuation theory;^"^ i.e., critical ex- 
ponents of thermodynamic and transport properties do 
not follow the standard theory, whereas in other cases 
the critical region is unexpectedly wide. In a number of 
compounds, this breakdown of the standard theory has 
been suggested in connection with the proximity of the 
first-order transition and the effects of inhomogeneities. 
For example, in weak itinerant ferromagnets ZrZn2,^ it 
has been proposed that the non-Fermi-liquid behavior 
is robust in a wide range of pressure. Similar behaviors 
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Fig. 1. (Color online) (a) Phase diagram with TCP under mag- 
netic fields. Continuous [First-order] phase transition line is rep- 
resented by thin (green) [thick (red)] curve, and the TCP is rep- 
resented by (yellow) circle, (b) Magnetization as a function of 
magnetic fields near the TCP. (1) At the first-order phase tran- 
sition, magnetization changes discontinuously. (2) At the TCP, 
magnetization changes non-analytically and its slope (FM sus- 
ceptibility xo = dM/dH) diverges. 




Fig. 2. (Color online) Magnetoresistance p{H) of YbRh2Si2 at 
2.3GPa reported in Ref. [17]. Discontinuous change of p{H) at 
low temperatures is the evidence of the first-order antiferromag- 
netic transition. 



are observed in MnSi^*^ and NiS2.^^ Furtliermore, near 
the first-order transition, a novel quantum pliase (ne- 
matic fluid) is found in Sr3Ru2 07^^ and unconventional 
superconductivity is found in UGe2.^'^ When continuous 
transition switches over to the first-order transition or 
phase separations, a tricritical point necessarily emerges 
as the boundary of these two. The purpose of this pa- 
per is certainly related to this motivation for elucidating 
physics under the proximity of first-order transitions and 
phase separations with the interplay with quantum fluc- 
tuations. 

A heavy- fermion compound YbRh2Si2 is a prototyp- 
ical example which does not follow the standard the- 
ory. It has been proposed that a diverging enhancement 
of uniform susceptibility xo occurs near the antiferro- 
magnetic (AFM) QCP.^'* The singularity of xo is esti- 
mated as Xo 2^^^ with ( ~ 0.6. Similar diverging en- 
hancements of the uniform susceptibility are observed in 
CeRu2Si2i^ (Xo oc T"? with C ~ 0.5) and in l3-YhAlBi^^ 
(Xo oc with ( ^ 0.3). One might speculate that this 
diverging enhancement could be caused by the hidden 
ferromagnetic (FM) QCP coexisting with the observed 
AFM QCP. However, the criticalities of these diverging 
uniform susceptibilities can not be explained by the con- 
ventional theory, because the critical exponent ^ must 
always be larger than one for the conventional FM QCP 
(see Table I). 

In our point of view, proximity to the first-order phase 
transitions is a key to understand the nature of these puz- 
zling quantum criticalities. Actually, in YbRh2Si2 and 
CeRu2Si2, evidences of the first-order AFM phase tran- 
sitions under magnetic fields are found^^'^* by tuning 
the pressure or substituting the chemical elements. In 
YbRh2Si2, at 2.3 GPa, it has been reported that resistiv- 
ity changes discontinuously as a function of the magnetic 
field at low temperatures (T < 0.5 K), while it changes 
continuously at high temperatures (T > 0.7 K)^'' as 
shown in Fig. 2. This indicates that the first-order tran- 
sition at low temperatures changes into the continuous 



one at higher temperatures through the tricritical point 
(TCP), as shown in Fig. 1 (a). In CeRu2Si2, although no 
clear magnetic order has been found in the stoichiomet- 
ric compound, AFM order appears in the Rh-substituted 
material Ce(Rui_a;Rh2.)2Si2 for x > 0.03. By apply- 
ing the magnetic fields, it has been observed that the 
first-order AFM transitions occur at low temperatures in 
Ce(Ruo.9Rho. 1)2812, while at high temperatures contin- 
uous AFM transitions occur. This experimental result 
indicates that the TCP also exists in Ce(Rui_j;Rha;)2Si2. 

In general, tricriticality necessarily induces an addi- 
tional divergence of uniform fiuctuations conjugate to the 
external field of the control parameter that drives the 
phase transitions. For example, under magnetic fields, 
the uniform susceptibility diverges at the AFM TCP. 

Here, we intuitively explain why the FM susceptibil- 
ity diverges at the AFM TCP under magnetic fields. As 
shown in Fig. 1(b), magnetization jumps at the first- 
order transition. By approaching the TCP, this jump be- 
comes smaller and vanishes at the TCP. Then, magneti- 
zation changes continuously but non-analytically at the 
TCP. This is the reason why the slope of the magneti- 
zation curve, namely, FM susceptibility diverges at the 
TCP. More precise discussions based on the (p^ theory 
are given in Refs. [19,20]. 

We note that the criticality of the classical TCP itself 
does not explain the unconventional quantum critical- 
ity observed in YbRh2Si2, CeRu2Si2, and /3-YbAlB4, be- 
cause the phase transitions are always either continuous 
or even absent, namely apparent phase transitions are 
not observed at ambient pressure and zero magnetic field. 
By taking YbRh2Si2 for example, we propose that the 
proximity effect of the quantum tricritical point (QTCP) 
is the possible origin of these unconventional quantum 
criticality. In YbRh2Si2, the first-order phase transition 
exists at high pressure, while the phase transition is 
always continuous at ambient pressure. From these ex- 
perimental results, by decreasing the pressure, we expect 
that the critical temperature of the TCP becomes zero 
and the QTCP appears at the critical pressure Pt [see 
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Fig. 3. (Color online) (a)- (c) Expected phase diagrams of 
YbRh2Si2 at various pressures. By decreasing pressure, it is ex- 
pected that the critical temperatures of the TCP become zero 
and the QTCP emerges at the critical pressure Pt. (a) Phase di- 
agram with the TCP [(yellow) circle] in T-H plane, where T (H) 
represents temperature (magnetic field). The TCP separates the 
continuous [thin (green) curve] and first-order [thick (red) curve] 
transition lines. In YbRh2Si2, similar phase diagram is proposed 
at 2.3GPai^ (see text), (b) Phase diagram with the QTCP. (c) 
Phase diagram of AFM phase with critical line [solid (green) 
curve] ending at the QCP [(yellow) diamond], (d) Phase diagram 
at zero temperature in P-H plane, where P represents pressure. 
The QTCP exists between the continuous and first-order transi- 
tion lines. 



Fig. 3(b)]. Because the QTCP is located very close to 
the ambient pressure, quantum tricriticality can be ob- 
served even at ambient pressure where the phase transi- 
tions are continuous. The criticality of the QTCP will 
be clarified in Sec. 3. In Sec. 4.1, we will show that 
the quantum tricriticality, which induces the divergence 
of uniform magnetic susceptibility, explains the uncon- 
ventional quantum criticality observed in YbRh2Si2. In 
Sec. 4.2 and 4.3, we will discuss that the proximity effect 
of the QTCP also explains the unconventional quantum 
criticalities observed in CeRu2Si2 and /3-YbAlB4. 

In this paper, we have clarified the criticality of the 
QTCP by extending the conventional spin fluctuation 
theory. In the conventional spin fluctuation theory, the 
origin of the non-Fermi liquid is ascribed to the coupling 
of the quasiparticle to the bosonic low-energy fluctua- 
tions of the order parameter. However, in the present 
quantum tricritical case, the quasiparticle couples not 
only to the bosonic order-parameter fluctuations but also 
to the uniform mode. Starting with this intuitive picture, 
we have shown that the serious modification of quantum 
critical phenomena arises from the equal and combined 
contribution of the two fluctuations which does not ex- 
ist in the conventional spin fluctuation theory. We note 
that the present spin fluctuation theory for the QTCP is 
applicable to the paramagnetic phase. 



A part of the spin fluctuation theory for the QTCP has 
already been briefly given in Ref. [21]. In this paper, we 
present the results of the quantum tricriticality in greater 
detail and discuss the singularities of the physical prop- 
erties near the QTCP, such as the uniform susceptibility 
Xo, the magnetization M, the specific heat 7, the nuclear 
relaxation time l/TiT, and the Hall coefficient i?H- We 
also give thorough comparisons with the experimental 
results. 

Before closing the Introduction, we briefly mention re- 
cent theory for the FIVI QTCP. The FJVI QTCP has been 
studied for itinerant helical ferromagnet MnSi'^^ by using 
the renormalization group theory. For nearly FIVI metal 
Sr3Ru207,^^ Green et al. have studied the QTCP by 
extending the SCR theory. However, we note that the 
singularity of Xq^ given by Green et al?^ as T^f^ is not 
correct, since they neglect the dependence of the bare 
second-order coefficient r^^ in their formalism. JVIoreover, 
these previous studies on the FJVI QTCP do not explain 
the unconventional coexistence of the FJVI and AFM fluc- 
tuations observed in YbRh2Si2.^^ 

The organization of this paper is as follows: In Sec. 2, 
we briefly review the conventional SCR theory. We 
mainly explain how the non-Fermi liquid behaviors ap- 
pear in the conventional SCR theory. In Sec. 3, we 
present the phenomenological SCR theory for the QTCP. 
From the present theory, we clarify the criticality of the 
QTCP; i.e., we obtain the critical exponents of the AFJVI 
susceptibility xq, the FJVI susceptibility xo, and the mag- 
netization curve. Section 4 describes comparisons of the 
present phenomenological SCR theory with the experi- 
mental results of YbRh2Si2, CeRu2Si2, and /3-YbAlB4. 
Section 5 is devoted to a summary and discussion. 

2. Conventional spin fluctuation theory for 
quantum critical point 

In this section, to make clear our starting point of the 
present study, we briefly review the conventional spin 
fluctuation theory for the quantum critical phenomena. 
The SCR theory proposed by JVIoriya is one of the stan- 
dard theory to describe the quantum critical phenom- 
ena. Originally, the SCR theory was proposed to ex- 
plain the weak and nearly FJVI or AFJVI metals. ^^"^^ Af- 
terwards, JVIoriya and Takimoto showed that the SCR 
theory can be applied to the quantum critical phenom- 
ena,'''^^ namely, they clarified how the spin fluctuations 
cause the non-Fermi liquid behaviors. We sketch the 
essence of the phenomenological SCR theory by following 
Refs. [6,7,28]. 

To understand the essence of the SCR theory, we start 
from a conventional Ginzburg-Landau- Wilson action for 
bosonic spin field ipq at the wave number q: 

SVPq] 



u 

iVo 



51 ■ "^-9') ^ (^9" ■ ^'3'-'?-'?"): (1) 

q,q',q" 

where u is a constant (we neglect the q dependence of 
the fourth coefficient u) and No is number of atoms. By 
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using this action, we obtain the free energy F as 

exp(-F/T)= [\{v^,eM-SW,]lT). (2) 

Because the contributions from the ordering wave num- 
ber Q are dominant for the conventional symmetry- 
breaking phase transitions, we approximate this free en- 
ergy as a function of AFM order parameter = (<Pq) 
up to the fourth order: 

f ~ i^o = Y^^' + + • • • ' (3) 

where fg is defined as rq = vq + 12uqIC, and spin fluc- 
tuation term /C is defined as 



(4) 



Here, we note that the bare second-order coefficient is 
renormahzcd by the spin fluctuation term IC. As wc will 
show later, this spin fluctuation term induces the non- 
trivial temperature dependence of the physical properties 
near the QCP. Hereafter, we only consider the paramag- 
netic state, i.e., = 0. 

By using the free energy in eq. (3), we obtain the or- 
dering susceptibility XQ as 



Xq 



dMr 



= fQ=rQ + UuqIC. (5) 

Mt^O 

According to the fluctuation-dissipation theorem,^^ the 
spin fluctuation term K, is described as 



q^Q 



2 f°° /I \ 



(6) 



where n(a;) = l/(e'^/-^ — 1). Near the QCP, we assume 
that x{Q + Qi ^) can be expanded with respect to q and 
iv as follows: 

X-i(Q + g,a;)~XQ'+ V-i^, (7) 

where 9 = l[e = 0] for the FM transitions {Q = 0) [AFM 
transitions (Q ^ 0)]. Here A and C arc constants. 
Prom cqs. (6) and (7), IC is described as 

IC = / dq duj{-+n{w)) 

Jo Jo 2 

d+e-i 



Ctjq" 



[q^iXQ^+Aq^W + iCco] 



+ 



— r dq r 

71" Jo Jo 
TT Jo Jo 



7-0-1 



[q'^iXQ^+Aq^W + iCur 
Cn{uj)ojq^+^-'^ 



[q^{XQ'+Aq^)]^ + {Cojy 



}C{0)+}C{T), 



(8) 



where d (T) represents the spatial dimensions (tempera- 
tures) and vq is the volume of the unit cell; Kj, is defined 



as Kd = Sd/{2TTY, where Sd = 2tt'^I'^ /T{d/2) and T is the 
gamma function; q^ is the cutoff wave number. Equations 
(5) and (8) constitute the self-consistent equation to be 
solved. 

For the brief notation, we introduce the variables as 



gs = 


{2d-K''-^/vo)-^, 


X = 




Ta = 


Aql 
2 ' 


To = 


A<&' 
27rC 






z' = 


w 


z 


27rr' 




2^' 


// 


2Ta ' 


t = 


T 

To 



Here, we note that the parameter Ta and To are so-called 
SCR parameters. Typical values of them in real materials 
are given in Refs. [6, 7]. By using these variables, we can 
describe /C(0) and IC{T) as 

Tod r , / r° x'^+^-^dx 
/C(0) = / dz' (9) 
-^A Jo Jo [{y + x^Jx^Y + z' 



IC{T) 



zdz 



2Tod 

~tZ Jo e^^--lJo [{y + x^)xyt]^ + z^- 

(10) 

Then, we obtain the singularity of /C(0) and /C(T). In 
three dimensions, because /C(0) can be expanded with 
respect to y regularly, we obtain the relation as 



/C(0) = Ko- Kiy. 



(11) 



Temperature dependence of IC{T) is obtained by the sim- 
ple scaling argument. To see the temperature depen- 
dence of /C(T) directly, we introduce the variable x' as 
X = x't^/'^^^ . By using this relation, we obtain 



K{T)=t^I{t), 



(12) 



Ta Jo e-^^-lh 



j-'d.i 



,<l+0-l 



[{y/t^ +x'^)x''']^ + z^ 
(13) 

In three dimensions, I{t) converges to a constant at zero 
temperature {t 0) if the condition 

hm y/t^/^+' ^ 

t— *oo 

is satisfied. 

In contrast to three dimensions, in two dimensions, 
I{t) has a logarithmic divergence and the analysis be- 
comes complicated. We do not show the complete anal- 
ysis here, because it is not the essential part of the SCR 
theory. To know the details of calculations for the two 
dimensional case, see Refs. [6,28]. Hereafter, we mainly 
consider the three dimensional case. 

From eqs. (5), (11), and (13), we obtain the relation 

y = yo + 12uQ{Ko-Kiy + t^I{t)), (14) 

where yo = rQ/2TA- This relation leads to 

d+e 

y = yo + yit^, (15) 

where yo = (yo + 12MQ_ft'o)/(l + i2uQKi) and yi = 
I{t)/{1 + UuqKi). The location of the QCP is given 
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Physical Properties 


3D AFM 


3D FM 


2D AFM 


2D FM 




^■3/2 


rp4/3 


-riog|iogr|/iogT 


-logT 


1 = C/T 


const. - 


-logT 


-logT 


J.-1/3 


1/TiT 


1/2 


XO 


XQ 


3/2 
Xo 


P 


^■3/2 


2^5/3 


T 


^4/3 



Table I. Critical exponents for the conventional QCP. For the AFM (FM) QCP, the ordering wave number Q is defined as Q 7^ 

{Q = 0). 



by yo = 0. Therefore, the singularity of the ordering sus- 
ceptibility near the QCP is obtained as 

Xq^ oc T=^/2 (3D AFM QCP), (16) 

Xo ^ oc T*/^ (3D FM QCP). (17) 

In two dimensions, by evaluating the singularities of /C(0) 
and /C(r) carefully, we obtain the singularities of order- 
ing susceptibility as 

Xq^ «-Tlog|logr|/logT (2D AFM QCP), (18) 

Xo^ « -T'logr (2D FM QCP). (19) 

By using the singularity of the ordering fluctuation, 
we can obtain the singularities of the specific heat,^° the 
nuclear relaxation time 1/TiT ,3i,32 ^^^^^ ^j^^ electronic 
resistivity p. '''33,34 Dgtails of the calculations are shown 
in Refs. [6,28]. Finally, we summarize the criticalities of 
the conventional QCP in Table. I. 

3. Spin fluctuation theory for quantum tricriti- 
cal point 

To clarify the criticality of the QTCP under magnetic 
fields, we extend the conventional Ginzburg-Landau- 
Wilson action in eq. (1), up to the sixth order of the 

bosonic field ip, as 

9 ^ g,Q',Q" 

X{(Pgn ■(pg,_g_g„) (20) 

+ 12 ('^91 ■ ^-92)(¥'93 • "P-Qi) 

^ifqs ■ </'q2+g4-9l -93-95) - H^O, (21) 

where H (Nq) is an external magnetic field (number of 

atoms); u{q, q' , q") and v are constants, while Vq depends 
on the magnetic field H . From eq. (21), the free energy 
F is obtained from 



where 



exp(-i^/T) = y[]P^,exp(-5[^,]/T). 



(22) 



Since the QTCP is expressed by fluctuations at both the 
AFM Bragg wave number q = Q and zero wave number 
q = 0, we approximate the free energy as a function of 
the order parameter M'l' = {(fq) and the uniform mag- 
netization M = (y^o)' 



+ -foM'^+UQM^ + vM^ -HM, (23) 



7 ^Q, ro, uo, and K. are defined as 
fQ{T, H) =rQ{H) + 12uq{1C + M^) 

uq{T,H) =(iQ + 15v(/C + m2), 
fo(T, H) =ro{H) + UuqIC + 90vlC^, 
uo{T,H) =uo + 15vK:, 



(24) 
(25) 
(26) 
(27) 

(28) 



95^0,Q 



Effects of spin fluctuations are included in K. fol- 
lowing the conventional SCR theory. We approximate 
■"(9, q, Q) [u{q, q, 0)] and the equivalent coefficients as q- 
independent values; u{q, q, Q) ~ uq [u{q, q, 0) ~ wq] for 
all q [for q ^ Q]- 

We eliminate M in eq. (23) by using the saddle-point 
condition for M, dFo/dM = 0, leading to the following 
relation between M and as 

M = ao + aiM^^ + a2M'''^ -\ , (29) 

where the expansion coefficients uq to 02 are determined 
by substituting eq. (29) into the saddle-point condition: 

fo(T, H)ao + 4{(o(T, H)4 + Gva^ ~ H = 0, (30) 

12aoUQiT,H)+aiR{T,H) = Q, (31) 

where R{T,H) = fo{T,H) + 12uo{T,H)al + SOva^ By 
using eq. (29), we obtain the free energy as 

Fo =lrQ{T, H)M^^ + u'q{T, H)M^^ + O(Mt'), (32) 

where u'q{T,H) = uq{T,H){1 + 6aoai). In eq. (32), 

continuous phase transitions occur at f q = when 
Uq{T,H) > 0, while the first-order phase transitions oc- 
cur when u'QiT,H) < 0}^''^° Therefore, the QTCP ap- 
pears when the conditions rQ{0, Ht) = and uq{0, Ht) = 
are both satisfied, where is the critical field at the 
QTCP. 

We now discuss the susceptibilities xq the AFM 
ordering vector Q and xo at g = in the disordered 
phase (M^ = 0, M = ao) by using eq. (30) and the free 
energy (32). From eq. (32), Xq^ is given as 

rQ{T,H). (33) 



Aft=o 



By differentiating eq. (30) with respect to the magnetic 
field H, we obtain Xo ^ as 

R{T,H) 



Xo 



1 j^Y^ = 

cxuq{T,H). 



1 - aodfo/dH - Aa^duo/dH 



(34) 



6 J. Phys. Soc. Jpn. 



Full Paper 



Online-Journal Subcommittee 



Here, we used eq. (31), which gives R{T, H) oc uq{T, H). 
Detailed derivation of eq. (34) is given in Appendix [see 
cq. (A-4)]. 

Here, we note that rQ{T, H) and uq{T, H) in the right 
hand sides of eqs. (33) and (34) are expressed by K. by 
using eqs. (24) and (25). Since K, is given from Imx by 
the fluctuation-dissipation theorem^^ as eq. (6), K can 
be expressed by using xq and Xo- Therefore, by us- 
ing eqs. (33) and (34), we can determine the singulari- 
ties of the susceptibilities (xg and Xo) Hclf-consistcntly. 
Once the singularities of xq and xo arc determined, we 
can determine the singularities of ao by using cq. (30). 
Hereafter, by using the above self-consistent equations, 
we clarify how the susceptibilities and the magnetiza- 
tion measured from the QTCP (xq, Xo, 5aQ = ao — aot 
with aot being the value at the QTCP) are scaled with 
6H = H-Ht and T near the QTCP. The results will be 
shown in eqs. (61)- (63) and Table II. 

As we have shown in the last section, in the SCR the- 
ory, the nontrivial temperature dependence of physical 
properties originates from the spin fluctuation term fC. 
Therefore, we need to clarify the scaling of K. by using 
the fluctuation-dissipation theorem by combining with 
expansions of Xo-fg(w) and XQ+q{^^) in terms of the wave 
number q and the frequency w near the QTCP. The 
fluctuation-dissipation theorem gives the relation 



(35) 

Hereafter, we mainly consider the three dimensional case. 

First, we consider the singularity of /C near the order- 
ing wave number. The ordering susceptibility XQ+q{^) 
is assumed to follow the conventional Ornstein-Zernike 
form, 

XQ+q{^)~^ - Xq^ + Aqq'^ - iCgUJ, (36) 

as in the conventional SCR formalism. From eqs. (35) 
and (36), in three dimensions, we evaluate /C near the 
ordering wave number as 



q^Q 



^ Jo Jo [{Xq 

^ Jo Vo "^[{XQ+AQq^Y + iCQU^Y] 



n{u))CQUjq^ 



= ICq{0)+ICq{T), 



(37) 



where /Cq(0) is the so-called zero point fluctuations and 
vq is the volume of the unit cell; is defined as = 
53/(277)3, where S3 = 2-^^/"^ /T{3/2) = An. For the brief 
notation, we introduce the variables as 



(67rVi^o)i 
27tCq' 



Tqa = 



Aqql 



Qbx, 



y{Q) 



xiQ)-' 



t{Q) 



T 
Tqo 



By using these variables, we obtain 

3Tqo /"^ . , f^" x^dx 



^q(O) = ^ / dz' 2 7 (38) 

^QA Jo Jo {y{Q) + x^y + z' 



6Jqo 
Tqa 



zdz r 

Jo e2--iyo 



X dx 



MQ)+x^)/tf+z-' 
(39) 



Here, we examine the temperature dependence of 
K-q{T) and /Cq(0). By scaling a; as a; = t^/'^x' , we ob- 
tain the singularity of K,q{T) as 



lCQ{T)=t"^lQ{t) 



(40) 



, 6T00 r zdz r 



x'^dx 



{y{Q)/t + x'y + z^ 



(41) 

Because self-consistent equations require the condition 
y{Q)/t ^ 1 for small t, Iqit) becomes constant at zero 
temperature. Therefore, near the QTCP, we obtain the 
singularity of /Cq (T) as 

K:q{T)(xT^/'^. (42) 

Because /Cq(0) can be expanded with respect to y{Q), 
we obtain the singularity of /Cq(0) as 

-1 



ICq{0) Kqo - Kc, 



(43) 



where Kqo and Kqi are constants. 

Next, we consider the singularity of K. near zero wave 
number. We note that the enhancement of the uni- 
form susceptibility is not caused by the conventional 
symmetry-breaking phase transition but is caused by the 
first-order AFM phase transition. Therefore, in contrast 
to the ordering susceptibility, the uniform part xa+q{'-^) 
does not follow the conventional Ornstein-Zernike form. 
Actually, according to the Ginzburg-Landau- Wilson the- 
ory, the scaling relation Xo ^(0) Xq^^^(O) holds near 
the TCP. As we will see, the self-consistency among 
eqs. (28), (30), (33), (34), and (35) requires that this 
relation still holds for g ^ 0. Therefore, we obtain the 
relation 

Xo+q{0)-' ocxQ+g(0)-V2 

ocixQ^+AQq^/'^iXo'+Aoq^^'. (44) 

Prom the conservation law, the ui dependence of 

Xo+g(w)"^ should be given as xo+<j(w)~^ ~ Xo+q(0)"^ - 
iCoOj/q. Finally, we obtain w and q expansions of 
Xo+q{u>)-'^ as 

Xo+qico)-' - (Xo"' + Aoq^/^ - tCoLo/q. (45) 
Prom eqs. (35) and (45), we evaluate K, near the zero 
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wave number as 



2 r°° /\ \ 



Physical Properties 


T dependence 


H dependence 


XQ 


rp-3/2 




X() 


rp-3/4 




5M 


rp3/4 





K3V0 



Jo "^Vo "^"^ [{Xo^ + Aoq^)q^ + {Coco)] 



H ^— ^ / dq I duj- j, _2 



n(cL')Coti'(?''' 



= /Co(0)+/Co(T), (46) 
For the brief notation, wc introduce the variables as 

x(o)-i 



y(o) 



T 

i(0) = 77^, 
J 00 



Too = qBToA/T^Co- 

By using these variables, we obtain 



J OA Jo Jo [y[Oy + x-^ix-^ + z' 



Jo Jo 



X dx 



67oo /"^ -2^-2 /"''^ x'^dx 

Tqa 



[(y(0)2 + a;2)/i]a;2 + z2 
(48) 



By scaling a; as a; = x'i}!'^, we obtain 

6T00 r zdz 



(49) 



To(t) 



x''^dx 



ToA Jo e2-^ - 1 Jo {y{0)yt + x'^)x'^ + ' 

(50) 



Because Io{t) becomes constant at zero temperatures, 
near the QTCP, we obtain the singularity of ICo{T) as 

/Co(T) oc T\ (51) 

Because ^^0(0) can be expanded with respect to y(0)2, 
we obtain the singularity of A^o (0) as 

ICoiO) ^ Koo - KoiXo^ (52) 

where JsToo and Kqi are constants. From eqs.(42), (43), 
(51), and (52), we obtain the singularity of (5/C = /C — /Cf 
measured from the critical value }Ct as 

SIC ~ -i^oiXo ' - KqiXq' + KotT^ + KqtT^/\ (53) 

Now the singularity of magnetization oq is obtained 
by solving eq. (30). Near the QTCP, eq. (30) can be 
approximated as 

A5al + BSao + C = 0, (54) 

with A = 12aot(5i;aQj + uo), B = 6fo + 12aQf6uo, and C = 
aotSfo+4:al^Sua - SH, where Svq = fo(T, H) - fo(0, Ht), 
and 6uo = uq{T,H) — uo{0,Ht). Since both B and C 
vanish at the QTCP, we obtain the asymptotic behavior 
of 6ao as 



Sao ^ {aoSH + aiSKf^ , 



(55) 



where ao and ai are constants. We give detailed deriva- 
tion of eq. (55) in Appendix [see eq. (A-1)]. 



Tabic II. Critical exponents for QTCP. Here, T and H represent 
temperature and magnetic field, respectively. 5M (SH) repre- 
sents the magnetization (magnetic field) measured from the crit- 
ical value. 



By defining 6rQ{T, H) = fgiT, H) - ^^(0, Ht), we ob- 
tain 

Xq' = 5fQ{T, H) =5rQ{H) + 9Qv{51C + 5~aof, (56) 

since both fQ(0, Ht) and uq{Q, Ht) are zero at the QTCP 
and terms linear in 5K and 5ao vanish. Here Srq and 5ao 
are defined as 

5rQ = TQiH) - TQ{Ht) - rgnSH, (57) 

Sclq = al - alt = 6aQ{5aQ + 2aoi). (58) 

From eqs. (53) and (55), 6IC turns out to be higher 
order of Sao near the QTCP. Then, from eqs. (34) and 
(56), the most dominant terms of Xg^ and Xo ^ are given 

as 

Xo ^ oc Sao, (59) 
Xq' ^ rQHSH + 360val,5al (60) 

which together with eqs. (53) and (55) lead to the 6H 
and T dependences as 

XQ'^fiQo5H + (3QiT^/', (61) 
Xo' ^ {PooSH + l3oiT^/y/\ (62) 

Sao^{a'„SH + a[T^/Y^^- (63) 

For simpler cases of (1) T 7^ and dH = (2) T = 
and SH ^ 0, readers are referred to Appendix for more 
detailed derivations of the self- consistent procedure. Sin- 
gularities of the uniform susceptibility xo, the ordering 
susceptibility XQ) and the magnetization Sao near the 
QTCP are summarized in Table II. Here, we note that 
these critical exponents are unique and any choice of the 
phenomenological parameters does not change the crit- 
ical exponents. These mean-field critical exponents ob- 
tained by the present SCR theory based on the classi- 
cal if^ theory are justified in three dimensions by the 
following reason: It is known that the upper critical di- 
mension dc is three ((ic=3) for the (p^ theory. For the 
quantum phase transitions, the effective dimension des 
is given as d^s = d + z, where z is the dynamical criti- 
cal exponent. For the AFM QTCP, since the dynamical 
exponent z is two, the effective dimension des is given 
as d(.g = d + z = 5, which is larger than the upper crit- 
ical dimension dc = 3. Therefore, the obtained critical 
exponents are correct in three dimensions. 

4. Comparison with experimental results 

We now examine whether the criticality of the QTCP 
is consistent with the experimental results for YbRh2Si2, 
CeRu2Si2, and /3-YbAlB4. In these three materials, the 
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convex temperature dependence of the inverse of the uni- 
form magnetic susceptibihty Xo^ (Xo ^ oc T'», < 1.0) is 
observed.^^'^®'^® We again note that this convex tem- 
perature dependence of Xo^ can not be explained by the 
conventional ferromagnetic quantum criticality, because 
the critical exponent ( is always larger than one. 

For YbRh2Si2, we will show that the present SCR the- 
ory reproduces quantitative behaviors of the experimen- 
tal uniform magnetic susceptibility and magnetization 
curve by choosing the reasonable set of the phenomeno- 
logical parameters. Furthermore, we show that the sin- 
gularity of the QTCP is consistent with the experimen- 
tal results of the specific heat, nuclear relaxation time 
1/TiT, and Hall coefficient. 

For CeRu2Si2 and /3-YbAlB4, we show that the sin- 
gularity of the diverging enhancement of the uniform 
susceptibility is consistent with the quantum tricriti- 
cality. We propose a possible location of the QTCP in 
CeRu2Si2. By finc-tmiing the Rh-substitution ratio and 
the magnetic fields, it is possible to determine the loca- 
tion of the QTCP. For /3-YbAlB4, we propose that NMR 
study is a suitable probe to verify the prediction of our 
QTCP scenario. 

4.1 YbRh2Si2 

4-1.1 Experimental results of YbRh2Si2 

At zero magnetic field {H = 0), it has been sug- 
gested that YbRh2Si2 exhibits an AFM transition at the 
Neel temperature Tn = 0.07 K.^^ Although neutron- 
scattering results are not available yet, anomalies of 
specific-heat and the uniform magnetic susceptibility in- 
dicate that the AFM transition actually occurs. Results 
of nuclear magnetic resonance (NMR)^^ and Mossbauer 
effect^® also indicate the existence of the AFM order. By 
applying the magnetic field along the ab plane, critical 
temperatures of AFM transitions become zero at the crit- 
ical magnetic field He - 0.06 T.^^'^^s p^j. Qe-substituted 
YbRh2(Sio.95Geo.o5)2, both Tn and He decrease to 0.02 
K and 0.027 T.^^ Because the AFM transitions still re- 
main continuous one down to the lowest temperature 
(~ 10m K), it has been proposed that a field-induced 
AFM QCP emerges. One might think that this clear 
QCP is a textbook example of the standard theory.'''"^ 
However, its non- Fermi-liquid properties do not follow 
the predictions of the standard theory and are under ex- 
tensive debates.^ 

The puzzling non-Fermi liquid behaviors observed in 
YbRh2Si2 are summarized as follows: At higher temper- 
atures (T > 0.3K), the Sommcrfeld coefficient of the 
specific heat 7 has logarithmic temperature dependence 
(7 oc — logT), while 7 is increased with power laws below 
0.3K."^^ This behavior is not consistent with the conven- 
tional theory, in which the 7 converges to a constant at 
low temperatures. Transport and optical data roughly 
show the resistivity measured from the residual resistiv- 
ity po, Ap = p — po, linearly scaled with T and fre- 
quency.^^' This singularity also contradicts the predic- 
tion of the standard theory (Ap cx T''/^). Moreover, it 
has been proposed that a large change in the Hall co- 
efficient Ru occurs near the QCP.^^ The most puzzling 
non-Fermi liquid behavior is a diverging enhancement 



of the uniform susceptibility xo near the AFM QTCP. 
The singularity of xo is roughly scaled by Xo and 
Xo \H — Hc\^'^ with ^ ~ ^' ~ 0.6^^ contradicting the 
standard expectation of saturation to a constant. NMR^^ 
and electron spin resonance (ESR)''° signals also indicate 
the enhancement of the FM susceptibility near the AFM 
QCP. In accordance with the diverging enhancement of 
the FM susceptibility, the magnetization curve has con- 
vex magnetic-field dependence. ^^'^^ It has been specu- 
lated that the origin of this enhancement could be the 
proximity to the FM QCP.^^ However, the critical expo- 
nent C ^ 0.6 can not be explained by the standard theory, 
because the critical exponent ( is always larger than one 
in the conventional FM QCP (see Table I). Under mag- 
netic fields H > He, it is proposed that two characteristic 
energy scales exist in YbRh2Si2;^^ one is the scale (Tlfl) 
for the establishment of the Landau Fermi liquid state, 
i.e., below Tlfl the resistivity has the Fermi- liquid form 
p = po + AT^, and the other one is the scale (T*) where 
dRii/dH, dp/dH, and xo have peaks. 

4-1-2 Choice of phenomenological parameters 

In this subsection, we explain how we choose the 
phenomenological parameters to solve the self-consistent 
equations in eqs. (28), (30), (33), (34), and (35) numeri- 
cally. Hereafter, for simplicity, we approximate the mag- 
netic field dependence of Sro{H) = ro{H) — ro{Ht) as 
6roiH)^roHSH. 

First, we clarify how many control parameters exist in 
the self-consistent equations. The four SCR parameters 
{ToA, Too, Tqa, and Tqo) are control parameters; the five 
parameters {v, rqn, toh, Ht, and aot) are also control pa- 
rameters, whereas once these parameters are fixed, the 
other parameters (ro, rq, uo, and uq) are determined 
from the conditions fQ{0,Ht) = 0, UQ{0,Ht) = and 
eqs. (30), (31). Therefore, the number of control param- 
eters is nine in the self-consistent equations. Below we 
show that these parameters are required to satisfy rather 
strict constraint from the physical reasons. 

Next, we explain how to employ a reasonable set of 
the control parameters. Because Ht and aot are the 
magnetic field and the magnetization at the QTCP, we 
can estimate these two parameters directly from experi- 
ments. In YbRh2Si2, since we expect that Ht and aot are 
slightly larger than those of the QCP at ambient pres- 
sure, we choose these two paramciter as Hi = 0.08T and 
aot = 0.2;Ub. According to the previous studies,^ it is 
suggested that the four SCR parameters (To^, Too, Tqa, 
and Tqo) have the values within the order of lO-lOOK. 
Therefore, we choose these parameters in this range; i.e., 
ToA = 20 K, Too = 10 K, Tqa=270 K, Tqo=5 K. In con- 
trast, for the other three non- primary parameters [rqu, 
roH , and v) , we do not find any constraint from the phys- 
ical requirement. Therefore, we have freely tuned these 
parameters to reproduce the experimental results quan- 
titatively {v = 21 K, roH = 5, and rQH = 180). How- 
ever, the critical exponents do not change even when we 
have chosen these parameters arbitrarily. In Appendix, 
we show the details of numerical calculations. The micro- 
scopic derivation of these phenomenological parameters 
is left for future studies. 
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4-1-3 Spin fluctuations and magnetization curve 

We now compare the numerical result of inverse of the 
uniform susceptibility Xo^ with the experimental result. 
As shown in Fig. 4, temperature dependence of Xo^ just 
above the QTCP is well consistent with experimental 
result. This consistency strongly supports the relevance 
of our proposal that the QTCP exists very close to the 
QCP in YbRh2Si2 as is shown in Fig. 3. Because the 
QCP in the experiment is located very close to the QTCP 
but slightly away from the QTCP, the nonzero offset of 
Xo^ exists in experiment. By approaching the QTCP, 
this nonzero offset becomes smaller and vanishes at the 
QTCP. 

In the last section, we have shown that the inverse of 
the uniform susceptibility Xo^ scales as T"-^^ in the low- 
temperature limit. As a first look, this scaling at asymp- 
totically low temperatures appears to be somewhat in- 
consistent with the experimental result {xo^ oc T"'^) re- 
ported in Ref. [14]. However, the numerical calculation 
shows that the asymptotic scaling deviates at finite tem- 
peratures and Xq^ looks nearly proportional to T*^-^ at 
higher temperatures (T > l.OK). This is consistent with 
the experimental result as shown in Fig. 4. 

Here, we discuss the competition between FM fiuctua- 
tions and AFM fluctuations near the QTCP. As shown in 
Fig. 4, at high temperatures (T > 0.6K), the magnitude 
of the FM fluctuations is larger than that of the AFM 
fluctuations (xo > Xq), while the magnitude of the AFM 
fluctuations is larger than that of the FM susceptibility 
at low temperatures (T < 0.6K). The origin of this ro- 
bustness of the FM susceptibility at high temperatures is 
the broad structures of xo+qi^) in the wavenumber space 
[see eq. (45)]. We note that this competition of the spin 
fluctuations near the QTCP is consistent with the exper- 
imental results of the NMR study on YbRh2Si2.^'* From 
the NMR study, it is proposed that the spin fluctuations 
are governed by the q = FM fluctuations away from 
the QCP, while near the QCP, the AFM fluctuations 
with a finite wave number far from q — develop sig- 
nificantly and governs the spin fluctuations. As we have 
shown here, this is indeed the tendency to be observed 
near QTCP. 

We also note that this competition of the spin fluctu- 
ations is a possible origin of the multiple energy scales 
observed in YbRh2Si2^'*''*^ as we mentioned in Sec. 4.1.1. 
In our QTCP scenario, T* is interpreted as the energy 
scale where the FM fluctuation xo begins to follow the 
Fermi-liquid form, while Tlfl is interpreted as the energy 
scale where the AFM fluctuation XQ begins to follow the 
Fermi-liquid form. Since both FM and AFM fluctuations 
diverge at the QTCP, two energy scales T* and Tlfl be- 
come zero at the QTCP. This result is consistent with 
the experimental result (see Fig. 2 B in Ref. 42). 

By using the same phenomenological parameters, we 
have also calculated the magnetization curve. As shown 
in Fig. 5, our result is well consistent with the exper- 
imental result in the interval of more than two orders 
of magnetic fields. We emphasize that the singularity of 
the magnetization curve {SM oc SH^^^, 6 = 2) can not 
be explained by the conventional FM quantum critical- 
ity, because the critical exponent S must be larger than 
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three for the conventional FM QCP. The present criti- 
cal exponent (5 = 2 is also completely different from that 
of the quantum critical end point, which belongs to the 
Ising universality class'^'^. For the Ising universality class, 
the critical exponent S is always larger than three. 




T(K) 

Fig. 4. (Color online)ExperimentaI uniform magnetic susceptibil- 
ity Xo^ for YbRh2(SiQ. 95660.05)2 at H = 0.03 T reported in 
Ref. [14] compared with numerical result of present SCR theory. 
Solid (red) [broken (blue)] curve represents the theoretical [ex- 
perimental] Xo^- Solid (green) curve represents the theoretical 
Xq^- The theoretical Xq^ and Xq^ are calculated just above the 
QTCP (H = Ht)- 




Fig. 5. (Color online)Experimental magnetization curve for 
YbRh2(Sio.95Geo.o5)2 at T = 0.09 K reported in Ref. [14] com- 
pared with the present theory. Solid (red) [broken (blue)] curve 
represents the theoretical [experimental] magnetization curve. 
SM (SH) represents the magnetization (magnetic field) measured 
from the critical value. We estimate the experimental critical 
magnetic field He (magnetization Mc) as 0.027 T (0.004 /^b)- 



4-1-4 Specific heat 

In the SCR theory, enhancement of x(g,w) is the ori- 
gin of the enhancement of effective mass (see Sec. 2 and 
Ref. [6,30]). Near the QTCP, because xili'^) has two 
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peaks around q = and q = Q, we can express the spe- 
cific heat as 

C _ 3K3V0 1-1^ 
T 



7 



71" JO 
■K 

IQ +70, 



dqq 



duj 



dqq^ 







duj 



c{uj)Ta+q 



(64) 



where c{lj) = cj^e^/^ /{T^e^/^ - 1)^), Tq+, = 
MXq'/Aq + q^)/CQ, and ro+, = AI^^XoV^o + 
q^y^^q/Co- By substituting the numerical result of XQ 
and Xo into this equation, we obtain the specific heat near 
the QTCP as shown in Fig. 6. We note that the contribu- 
tion from zero wave number (70) is comparable to that 
of the ordering wave number (7q). This indicates that 
the quantum tricriticality induces larger enhancement of 
the effective mass than that of the conventional quantum 
criticality. This might be the reason why YbRh2Si2 has 
larger effective mass'^^ (7 ^ 1.5J mol^^K^^) than those 
of other typical heavy-fermion compounds (7 < 1.0 J 
mol-iK-2). 

Here, we discuss the obtained singularity of the Som- 
merfeld coefficient of the specific heat 7. At high temper- 
atures {T > 1.0 K), the singularity (7 oc — logT) and the 
amplitude are both consistent with those of the experi- 
mental result. However, at low temperatures (T < 1.0 
K), within this SCR theory, the singular temperature 
dependence of 7 is the same as that of the conventional 
AFM QCP (7 oc const. — T^/^), while experimentally, 
power-law-like behavior is observed for T < 0.3 K.'^^ In 
general, as long as we consider the spin fluctuations, we 
can not obtain the power-law divergence of the specific 
heat in three dimensions. 

Although the clarification of the origin of this discrep- 
ancy is left for future studies, we point out the two pos- 
sible origins of this discrepancy. One is the fact that the 
Neel temperature is actually nonzero in the experimen- 
tal results. At zero magnetic field, specific heat has a 
sharp peak around the Neel temperature.'^^ Remnants of 
this peak may be the origin of the power-law-like behav- 
ior of the specific heat at H = 0.06T. The other pos- 
sible origin is effects of valence fiuctuations,^^''*^ which 
is not considered in our spin fluctuation theory. In gen- 
eral, the quantum tricriticality, namely proximity to the 
first-order AFM transition, induces the divergence of the 
valence fiuctuations by the following reasons: Discontinu- 
ous change in the occupations of / electrons occurs at the 
first-order phase transition, and at the QTCP, it is ex- 
pected that the valence of / electrons changes singularly 
and the valence susceptibility diverges. We note that, in 
contrast to the pure valence transition, this divergence 
of the valence fluctuations is caused by the proximity to 
the first-order AFM transition. Therefore, in this case, 
the criticality of the valence fluctuations is governed by 
the quantum tricriticality. In other words, the Ising crit- 
icality expected from the simple valence quantum insta- 
bility should not show up in this case, similarly to our 
statement that the uniform susceptibility does not follow 
the ferromagnetic criticality. However, it is expected 



that contribution to the specific heat from the valence 
fluctuations, which is not considered in our spin fluctu- 
ation theory, explains the anomalous enhancement of 7 
observed in experiment. 




r(K) 



Fig. 6. (Color online)Numerical result for the Sommerfeld con- 
stant of the specific heat just above the QTCP. Here, 7q(7o) 
represents the contributions from the ordering (zero) wave num- 
ber. 



.1.5 NMR 
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Fig. 7. (Color online)Singularity of the 1/TiT just above the 
QTCP. We note that l/TiT looks proportional to T'^/^ at 
the intermediate temperatures (10~^K < T < 10~^K). Near 
the QCP in YbRh2Si2, 1/T\T appears to be proportional to 
2^-1/2 24 'pjjjg behavior is consistent with the quantum tricriti- 
cality. 



In this section, we consider the singularity of the nu- 
clear relaxation time 1/TiT near the QTCP. As we have 
mentioned in Sec. 2, the singularity of 1/TiT is scaled 
with the spin fluctuations.'^'^' By using the conventional 
relation, we obtain the singularity of 1/TiT near the 
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QTCP as 
TiT No '^'^^0 uj 

^Jo ^0+q Jo ^Q+q'^ 

^ -DologT + DqT-''/^ (65) 

where 7jv is a gyromagnetic ratio and Aq is the q compo- 
nent of hyper-fine couphngs; Do and Dq are constants. 

We note that the contributions fi-om FM susceptibihty 
induce the unconventional logarithmic temperature de- 
pendence of l/TxT. Because it is difficult to determine 
the value of Aq and Ao quantitatively, we do not calcu- 
late to show quantitative behavior of the 1/TiT. Alter- 
natively, in Fig. 7 (a), we show a qualitative behavior of 
l/TiT by simply taking Do = Dq = 1. Because the sub- 
dominant term — logT exists, l/TiT does not follow the 
dominant temperature dependence T~^^'^ at high tem- 
peratures (T > 10~^ K) but seems to be proportional to 
T~^^'^ at the intermediate temperatures. This behavior is 
consistent with the unconventional temperature depen- 
dence (l/TiT oc observed in YbRhaSia.^'' 

4-. 1.6 Hall coefficient 

In this section, we show that the singularity of the Hall 
coefficient can be explained by the quantum tricriticality. 
According to the Coleman's simple argument in Ref. [46], 
the Hall coefficient i?H is proportional to the square of 
the order parameter (i?H oc iW^^).® Near the QTCP, 
the order parameter is proportional to \g — gc\^^^ (see 
Refs. [19,20]), where g {gc) is the (critical value of) con- 
trol parameter. Therefore, Hall coefficient is scaled with 
\g — gc\'^^'^ near the QTCP. This result indicates that the 
Hall coefficient changes non-analytically near the QTCP. 
Moreover, if the QCP in YbRh2Si2 is located slightly 
away from the QTCP but on the side of weak first-order 
phase transitions, the Hall coeSicient must change dis- 
continuously with a jump. These behaviors are consistent 
with the experimental results. 

We note that the above argument is qualitative. To 
clarify the change in the Hall coefficient quantitatively, 
it is necessary to consider the realistic band structures 
of YbRh2Si2. By calculating the band structures of 
YbRh2Si2, Norman^^ has pointed out that small changes 
of the / electron occupation are sufficient to reproduce 
the experimental result. In general, proximity to the first- 
order transition indeed induces such changes of the / 
electron occupation near the QTCP. Although it is in- 
triguing to perform detailed calculations of the Hall co- 
efficient based on the microscopic band structures com- 
bined with the quantum tricriticality, such treatment is 
beyond the scope of the present paper. 

12 CeRu2Si2 

CeRu2Si2 is a canonical heavy-fermion system, which 
has no apparent magnetic order and shows a large Som- 
merfeld coefficient 7 ~ 360mJ/molK^.'**'''^ The result of 
the neutron scattering shows that AFM spin correlation 



develops below 60K,^° and detailed inelastic neutron- 
scattering study shows that the spin fluctuations can be 
explained by the conventional SCR theory. '"^^ By substi- 
tuting Ru (Ce) with Rh (La) slightly,^^'^^ AFM long- 
range order emerges. From these experimental results, it 
is recognized that CeRu2Si2 is located very close to the 
AFM QCP. We note that, although no clear evidence of 
magnetic long-range order was observed by the measure- 
ment of the bulk properties, ultrasmall ordered moment 
(- IO^Vb/Cc) was detected by ^SR below O.IK.^^ 

In CeRu2Si2, specific heat, resistivity, and uniform 
magnetic susceptibility do not show non-Fermi liquid 
behaviors down to 20mK. However, below 20mK, non- 
Fermi liquid behavior is observed. Actually, as shown in 
Fig. 8, Takahashi et al. have found a diverging enhance- 
ment of the uniform magnetic susceptibility at tempera- 
tures of micro Kelvin order. In the same temperature 
range, Yoshida et al. have found that thermal expan- 
sion and magnetostriction also show non- Fermi liquid be- 
haviors.^* By fitting the uniform magnetic susceptibility 
with the function Xo^ = aT'^, we estimate the critical 
exponent C a-s ^^0.53. This value of C is close to that 
of the QTCP, and is not explained by the conventional 
quantum criticality. Therefore, it is plausible that this 
diverging enhancement is caused by the proximity effect 
of the QTCP. 




I ^ ^ ^ ^ 1 

1 2 3 4 5 

r(mK) 

Fig. 8. (Color online)Experimental Xq^ ior CeRu2Si2 at low tem- 
peratures reported in Ref. [15]. The solid line shows the result 
of a least-square fitting by assuming the function = aT^, 

where a is a constant. We estimate the critical exponent ( as 
C ~ 0.53. 

We now discuss where the QTCP is lo- 
cated in CeRu2Si2. In Rh-substituted compound 
Ce(Rui_a;Rh2;)2Si2, it is known that the AFM magnetic 
order occurs for x > 0.03. At x—0.1, it was shown that 
the slope of the magnetization curve becomes steeper 
by lowering temperatures (see Fig. 9 in Ref. [18]). This 
experimental result indicates that the continuous AFM 
phase transitions changes into the first-order ones at low 
temperatures, and between them TCP exists at finite 
temperatures, where the slope of the magnetization 
is supposed to diverge. We show the expected phase 
diagram for x — 0.1 in Fig. 9(a). We note that this 
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(SL)X>Xt (h)x = xt 




H=Ht 



Fig. 9. (Color online)(a) Phase diagram of Ce(Rui_3;Rha;)2Si2 
for X > xt. Here, xt is the critical substituting ratio where the 
critical temperatures of the TCP becomes zero, (b) Phase dia- 
gram of Ce(Rui— xR.h2:)2Si2 at a; = xt. At this substituting ratio, 
the QTCP appears a.t H = Ht- (c) Expected phase diagram of 
Ce(Rui_a;Rha;)2Si2 at zero temperatures. The QTCP appears 
at (x,H) = {xt,Ht). 



first-order magnetic phase transition is not related with 
the metamagnetic transitions at Hyi ^ 8T observed 
in pure CeRu2Si2. The remnant of this metamagnetic 
phase transition still survives at x = 0.1 around H ~ 6T. 
By decreasing the substitution ratio x, it is expected 
that the critical temperatures of the TCP become 
zero at the critical substitution ratio xt and QTCP 
emerges as is shown in Fig. 9(b). We show the expected 
ground-state phase diagram of Ce(Rui_2:Rlia;)2Si2 in 
Fig. 9(c). Similarly to YbRh2Si2, because the QTCP 
is very close to x — 0.0 and H — 0.0, the diverging 
enhancement of xo is observed in pure CeRu2Si2. 

Here, we comment on the singularities of other phys- 
ical properties (specific heat, nuclear magnetic relax- 
ation time, and Hall coefficient), which are discussed for 
YbRh2Si2. In these physical properties, it is expected 
that the same singularities as those of YbRh2Si2 are ob- 
served in CeRu2Si2, because the criticality of the QTCP 
does not depend on the details of materials. 

Compared with YbRh2Si2, the occurrence of the di- 
verging enhancement of xo is suppressed to very low 
temperatures (T < 20mK) in CeRu2Si2. Although the 
criticality of the QTCP accounts for the singularity of 
Xo qualitatively, it does not explain the origin of this 
suppression. Clarifying the origin of this suppression is 
left for future studies. 

To summarize, in CeRu2Si2, QTCP is expected to ex- 
ist very close to the ambient pressure and the zero mag- 
netic field, and quantum tricriticality is the origin of the 



anomalous diverging enhancement of the uniform mag- 
netic susceptibility. It is highly desirable to determine 
the location of the QTCP precisely by fine-tuning the 
Rh-substitution ratio x and the magnetic field H . 

4.3 P-YbAlBi 

/3-YbAlB4 is a newly discovered heavy-fermion com- 
pound, which has a large Sommerfeld coefficient 7 
300mJ/molK^.^^ In this material, up to now, neither ap- 
parent QCP nor magnetic order have been found at am- 
bient pressure. However, the superconducting transition 
with Tc = 80mK is found.^^.se 

In /3-YbAlB4, the non- Fermi liquid behaviors are ob- 
served at ambient pressure. The striking feature of (3- 
YbAlB4 is the divergence of the uniform magnetic sus- 
ceptibility (xo oc r^^, C = 1/3).^^ Because the critical 
exponent ^ is smaller than one, we can exclude the pos- 
sibility that this divergence is caused by the FM QCP. 
However, this critical exponent = 1/3 is substantially 
smaller than that of the QTCP. If this value 1/3 is 
correct, a possible origin of this discrepancy is the low 
dimensionality. From the crystal structures, it is ex- 
pected that this material has one dimensional anisotropy. 
This low dimensionality may modify the criticality of the 
QTCP at high temperatures (T > O.IK). However, the 
criticality is expected to follow that of the three dimen- 
sionality at sufficiently low temperatures {T ^ O.IK). 
We expect such a crossover to occur at lower tempera- 
tures below O.IK. 

Here, we explain why the low dimensionality makes the 
critical exponent C small. From eq. (15), critical exponent 
of the order-parameter fluctuations becomes smaller by 
lowering the dimensions (see also Table I). Therefore the 
critical exponent C, which is scaled with half of the crit- 
ical exponents of the order-parameter fluctuations, also 
becomes smaller by lowering the dimensions. 

If the AFM QTCP exists around the ambient pressure 
and the zero magnetic fleld in /9-YbAlB4, coexistence of 
the enhanced AFM and FM fluctuations are expected 
to be observed (see Sec. 4.1.3). By performing the NMR 
measurement, it is possible to detect such coexistence of 
spin fluctuations as is observed in YbRh2Si2.^^ 

5. Summary and Discussion 

In this paper, quantum tricriticality has been stud- 
ied by extending the conventional SCR theory. Near the 
AFM QTCP, not only the AFM fluctuations but also the 
FM fluctuations show diverging enhancement. By con- 
sidering the combined effects of these two different di- 
verging fluctuations, we have shown that unconventional 
non-Fermi-liquid behaviors appear around the QTCP. 
We have proposed that the quantum tricriticality ex- 
plains the unconventional quantum criticality observed 
in YbRh2Si2, CeRu2Si2, and /3-YbAlB4. 

In YbRh2Si2, as we have explained in Sec. 4.1, avail- 
able experimental results strongly support the existence 
of the QTCP under pressure. Actually, we have shown 
that the quantum tricriticality is consistent with the puz- 
zling quantum criticality observed in YbRh2Si2. 

More concretely, near the AFM QTCP, we have clar- 
ified singularity of the uniform magnetic susceptibil- 
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ity as xo k T""^/"* at low temperatures. We have also 
clarified the singularity of the magnetization curve as 
5M oc SH^f^. It is noteworthy that these critical ex- 
ponents are completely different from the conventional 
quantum criticality (see Table I), and are consistent with 
the experimental results of YbRh2Si2. Furthermore, for 
the magnetic susceptibility and the magnetization curve, 
by solving the self-consistent equations numerically, we 
have shown that the quantum tricriticality is consistent 
with the experimental results of YbRli2Si2 not only qual- 
itatively but also quantitatively. This includes a crossover 
from Xo «: T~^/^ to xo ~ T~°-^ at elevated tempera- 
tures. 

On the Sommerfeld constant of the specific heat 7, 
we have shown that the singularity (7 oc — log T) and 
the amplitude arc both consistent with those of the ex- 
perimental result'^^ at high temperatures (T > l.OK). 
Especially, near the QTCP, we have shown that the con- 
tribution from the uniform magnetic susceptibility plays 
an important role in making the heavy quasiparticles. 
However, at low temperatures (T < l.OK), within our 
tlicjory, 7 converges to a constant (7 oc const. — T^^"^), 
while experimentally, a power-law-like behavior is ob- 
served for T < 0.3K.^^ This discrepancy may be solved 
by considering either the fact that the Neel temperature 
is actually nonzero or effects of valence fluctuations.^''' 
In our point of view, contribution from the valence fluc- 
tuations is the most promising candidate to explain the 
anomalous enhancement of effective mass. It is an in- 
triguing challenge to extend the present spin fluctuation 
theory to treat the interplay between spin and charge 
fluctuations such as valence fluctuations. 

The quantum tricriticality also induces unconventional 
behavior of the nuclear relaxation time, namely 1/TiT. 
Near the QTCP, we have shown that the 1/TiT is pro- 
portional to -Do log T + DqT-^/^, where Dq and Dq 
are constants. Due to the existence of a logarithmic di- 
vergence, temperature dependence of the 1/TiT seems 
to be weaker than T~^/^ at high temperatures. Thus, at 
intermediate temperatures, we have shown that l/TiT 
seems to be proportional to T~^l'^ . This behavior is con- 
sistent with the experimental results of YbRh2Si2.^^ 

Wc have also shown that the large change in the Hall 
coefficient^^ can be well accounted for by the quantum 
tricriticality. In this paper, we only consider the quali- 
tative aspect of the Hall coefficient. It is left for future 
study to clarify the quantitative changes of the Hall co- 
efficient near the QTCP. 

Here, wc comment on singularities of the resistivity. In 
the quantum tricritical scenario, if we consider the sin- 
gularity of the relaxation time only, the singular temper- 
ature dependence of the resistivity is the same as that 
of the conventional three dimensional AFM QCP, i.e., 
p oc T^/^. This is inconsistent with the linear-like tem- 
perature dependence of the resistivity observed in experi- 
ment.'^^ However, according to the simple Drude picture, 
the resistivity p is proportional not only to the relaxation 
time T but also to the carrier density n, i.e., p oc l/nr. 
Near the QTCP, it is expected that the carrier density 
also has the singularity. If its singularity is the same as 
that of the magnetization, we obtain the singularity of 
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the resistivity as p oc T'^/^/ (ho + niT^/^), where no and 
ni are constants. This singularity of the resistivity may 
explain the linear-like temperature dependence of the re- 
sistivity observed in the experiments. To examine the va- 
lidity of the above simple arguments, further studies are 
desirable. 

As we have explained above, many of unconventional 
behaviors of YbRh2Si2 can be accounted by the quantum 
tricriticality. In fact, we do not find any experimental re- 
sults of YbRh2Si2 that explicitly contradict our theory. 
We emphasize that no other theory in the literature is 
able to explain the anomalous quantum criticality ob- 
served in YbRh2Si2 such as the coexistence of the en- 
hanced AFM and FM fluctuations, and the criticality of 
the uniform susceptibility. Although these results ensure 
the relevance of our theory, further experimental study 
to determine the precise location of the QTCP by tuning 
the pressure and the magnetic field will be a crucial test 
of our theory. 

Here, we comment on three different scenarios pro- 
posed for the QCP in YbRh2Si2; local quantum crit- 
icality,^^' reconstruction of the Fermi surface, and 
fermion condensate. 

In the local quantum critical scenario, Coleman et al. 
claim that a breakdown of a composite heavy fermion 
(namely, all / electrons decouple from the Fermi surface) 
occurs at the AFM QCP. In a simple interpretation, their 
scenario indicates that no heavy electron exists in the or- 
dered phase any more. Thus, their scenario seems to be 
inconsistent with the large Sommerfeld coefficient of the 
specific heat observed in the ordered phase. Although 
they claim that the dynamical Kondo correlations in the 
ordered phase can account for the heavy electrons,^' 
to the authors' knowledge, there are no quantitative cal- 
culations to be compared with the experimental results. 
It is desirable to examine whether a large Sommerfeld 
coefficient of the specific heat 7 observed even in the or- 
dered phase can be quantitatively reproduced in the local 
quantum critical scenario. 

We now examine the expcrinicintal evidence of the lo- 
cal quantum criticality. It has been proposed that a large 
change in the Hall coefficient in YbRh2Si2^^ is the evi- 
denc;e of the local quantum criticality. However, as we ex- 
amined in Sec. 4.1.6, quantum tricriticality naturally ex- 
plains such a large change in the Hall coefficient qualita- 
tively. Therefore, the large change in the Hall coefficient 
is not a conclusive evidence of their scenario. The singu- 
larity of the nuclear magnetic relaxation time \/TiT is 
given as 1/TiT oc in their local quantum critical sce- 
nario (see eq. (4) in Ref. [57]). We emphasize that this 
behavior is not consistent with the experimental result 
of the NMR measurement.^* The local quantum criti- 
cal scenario also claims that the susceptibility follows 
the Lo/T scaling. Although this behavior is observed in 
CcCus.gAuo.i^^ , it is not clear whether this behavior can 
be observed in YbRh2Si2 because the neutron-scattering 
results are not available yet. In contrast to their sce- 
nario, our theory predicts the conventional tj/T^/^ scal- 
ing.^ Further experiments are desirable to examine which 
behavior is observed in YbRh2Si2. 

Here, we note that the unconventional quantum criti- 
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cality observed in CeCu5.gAuo.i is qualitatively different 
from that of the YbRh2Si2. In particular, no diverging 
enhancement of the uniform magnetic susceptibility is 
observed and no clear evidence of the first-order transi- 
tion is found. Therefore, quantum tricriticality is not an 
origin of the unconventional quantum criticality observed 
in CeCu5.gAuo.i. Although the local quantiun critical- 
ity can explain the inelastic neutron-scattering result^^ 
and large changes in the Hall coefficient,^^ it can not ex- 
plain all the c;xperimental results. For example, singular- 
ity of the nuclear magnetic relaxation time 1/TiT does 
not simply follow the prediction of the local quantum 
criticality.^"^' Furthermore, by changing the Au substi- 
tution ratio slightly, it has shown that inelastic neutron- 
scattering result is described better by the conventional 
uj/T^/"^ scaling than by the lj/T scaling.^'"'' These exper- 
imental results indicates that the local quantum crit- 
icality is insufficient for explaining the unconventional 
quantum criticality of CeCus.gAuo.i- Further studies are 
needed to clarify the origin of the unconventional quan- 
tum criticality observed in CeCus.gAuo.i. 

Next, we comment on the Fermi-surface reconstruction 
scenario. By using the variational Monte Carlo method, 
Watanabe and Ogata^^ claim that the first-order AFM 
transition accompanied by the changes in a Fermi-surface 
topology occur in the two-dimensional Kondo lattice 
model. Their scenario is different from the local quan- 
tum criticality in the sense that the Kondo screening 
still remains even in the AFM phase. They also claim 
that this first-order phase transition is the possible ori- 
gin of the large change in the Hall coefficient observed in 
YbRh2Si2. Their scenario has similarity with our quan- 
tum tricritical scenario on the point that the proximity to 
the first-order phase transition induces the large change 
in the Hall coefficient. However, we note that the exper- 
imental large change in the Hall coefficient is observed 
near the continuous phase transition. Therefore, it is not 
clear whether the nature of this first-order phase tran- 
sition can explain the experimental large change in the 
Hall coefficient. Further study on the interplay between 
the changes in Fermi-surface topology, i.e., Lifshitz tran- 
sitions^^"™ and the quantum tricriticality is a challeng- 
ing issue for the future. 

In the fermion-condensate scenario, by assuming the 
divergence of the effective mass, it is shown that both 
the FM susceptibility xa the Sommerfcld coefficient 
of the specific heat 7 diverge with the singularity T~^^'^. 
They proposed that this singularity is consistent with 
the experimental results of YbRh2Si2 and CeRu2Si2. Al- 
though their theory has succeeded in explaining the ex- 
perimental results, it is not clear why the effective mass 
diverges at the AFM QCP. Furthermore, we point out 
that their theory does not explain why the FM and AFM 
fluctuations coexist as observed in YbRh2Si2. 

In Sec. 4.2, we have pointed out that the quantum 
tricriticality can be observed in CeRu2Si2. In CeRu2Si2, 
the diverging enhancement of the FM susceptibility xo 
is observed at very low temperatures.^''' We have shown 
that the singularity of xo is qualitatively consistent with 
the criticality of the QTCP. This result indicates that the 
QTCP is located very close to pure CeRu2Si2. For the 



Rh-substitutcd compound Ce(Rui_2:Rh2:)2Si2, we pre- 
dict that the AFM TCP exists under magnetic field at 
X = 0.1. By decreasing the substitution ratio x, it is ex- 
pected that the critical temperatures of TCP becomes 
zero and the QTCP appears at the critical substitution 
ratio Xt (see Fig. 9). CeRu2Si2 appears to be a suitable 
material to search for the QTCP, because large single 
crystals are available and precise measurements are pos- 
sible. 

In Sec. 4.3, we have proposed that the quantum tri- 
criticality is the possible origin of the diverging enhance- 
ment of Xo observed in /3-YbAlB4. Although the convex 
temperature dependence of Xo^"'^ is consistent with the 
criticality of the QTCP, the critical exponent C = 1/3 is 
much smaller than that of the QTCP. We have pointed 
out that a possible origin of this discrepancy is the low 
dimensionality of /3-YbAlB4. It is left for future studies 
to clarify how the low dimensionality affects the critical- 
ity of the QTCP at finite temperatures. 

We now discuss the singularities of the FM suscep- 
tibility around the conventional AFM QCP. Based on 
the conventional SCR theory, Hatatani, Narikiyo, and 
Miyake clarified the singularity of the FM susceptibility 
around the AFM QCP at zero magnetic field with the dy- 
namical exponent z = 2 as well as 2; = SJ^ They showed 
that the singularity of xo is given as xo ~ a — 6T^/^ 
(Xo ~ a - &rV3) for ^ = 2 (z = 3) AFM QCP, where a 
and b are constants. They proposed that such singular- 
ity is consistent with the experimental result of CerNia.''^ 
By using the renormalization-group theory, Fischer and 
Rosch proposed that the FM susceptibility has the sim- 
ilar singularity around the AFM QCP under magnetic 
fields. We note that these theories are only applicable 
to the continuous quantum phase transitions, and do not 
reproduce the diverging enhancement of the FM suscep- 
tibflity observed in YbRh2Si2, CeRu2Si2, and /?-YbAlB4. 

Finally, we discuss broader implications of the quan- 
tum tricriticality for the other fields of condensed matter 
physics. We have applied the theory of quantum tricriti- 
cality specifically to typical / electron systems with large 
effective masses. However, it should be noted that the 
present theory offers a general framework which may be 
applied to systems other than the heavy fermion com- 
pounds. The quantum tricriticality is significant in sev- 
eral different fields of condensed matter physics and will 
attract much interests from experimentalists as well as 
from theorists, because the TCP is known to play im- 
portant roles in various fields of physics, for instance, in 
the problem for the mixture of "'He and ^Hc^'* as is stud- 
ied by using Blume-Emery-Griffiths model. Recently, 
it has also been proposed that the TCP plays an impor- 
tant role on ultracold atomic Fermi gases.™ Therefore, 
it is a highly fundamental issue how the tricritical phe- 
nomena are modified when quantum fluctuations drives 
the critical point to zero temperature. It is desirable to 
explore the novel phenomena such as unconventional su- 
perconductivity induced by the quantum tricriticality. 
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Appendix: Details of the numerical calculations 
for the self-consistent equations 

In this appendix, we show the details of the numeri- 
cal calculations for the self-consistent equations for the 
QTCP given in Sec. 3. For simplicity, we only consider 
the solutions for two cases; (1) T 7^ and SH = 0, and 
(2) T = and 5H ^ 0. 

(1) T ^ and dH = 

First, to determine the singularities of the physical 
properties, we will obtain the five self-consistent equa- 
tions given in (A-l)-(A-5). By solving eq. (54), we obtain 
the singularity of 6ao as 



defined as 



5ao 



2At. 



(A-l) 



where wc approximate A — 12aQt{5vaQ^ + uq) a.s A ^ 
At = 12{(otaot-|-60uaQj, and B (C) is defined as i? = Sro+ 
12a§j(5uo (C = aot^fo + ^ciot^uo — dH). From eq. (56), we 
obtain the singularity of Xg^ as 

Xq' = SfQ{T, H) = SrQ{H) + 90v{SIC + Saof. (A-2) 

We also obtain the singularity of SK. as 

SIC = -KoiXo^ - KqiXq' + /Co(T) + )Cq{T). (A-S) 

By differentiating eq. (54) with respect to the magnetic 
field H, we obtain the singularity of xo as 



Xo = - 



9B s:„ I dC 

25aQAt + B 



(A-4) 



where 

Ml = -[{l%Qv8K, + At/aot + IZOwa^J 

5ao + At + 180t;aot]/(2^ao^t + B), 

M2 = -{ronSao + aotroH - l)/{25aoAt + B). 

Here, wc note that xo depends on dSIC/dH explicitly. To 
obtain the explicit form of the dSIC/dH, we differentiate 
SK. with respect H; 

= -i^Ol—^^TT- H 



dH 



Kqi 



dH ■ dH dxo' 
9Xq' dxQ'dlCQ{T) 

dH ^ dH dXg' 

dJCo{T) 



-Bo{K 



01 



dXo- 



)-Bq{Kq, 



dx7' 



Q 



(A-5) 

where we approximate dx^^/dH and dxq^/dH as their 
values in the low-temperature limit, i.e., Bq and Bq are 



Bq = lim 



dXo'- 



Bn = lim 



T^o dH 



T^o dH 



6H=0 



Bq and Bq can be determined from eqs. (A-l)-(A-4). By 
solving eqs. (A-l)-(A-5), we can determine the singularity 
of 

Xo ^> Xq^! ^oiO! <5^! and dSlC/dH. Concrete procedure 
is shown in Fig. A-l. 
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Fig. A-l. Schematic diagram for calculating self-consistent equa- 
tions given in eqs. (A-l)-(A-5). Here, e is the criterion for the 
convergence. 



(2) T = and SH ^ 

At zero temperature, from eq. (A-3), we can easily ex- 
press Xo as a function of 6)C because both /Co(0) and 
/Cq (0) are zero; 



Xo — I 7^ ) 

Aoi 



(A-6) 



We note that Xo ^ only depends on SIC and SH, because 
Xq^ is the function of SK. and SH (sec eqs. (A-l) and 
(A-2)). Therefore, by using eqs. (A-4) and (A-6), we ob- 
tain the differential equation with respect SIC as 

9SIC , , X ^ / , ^ 
(Xo - M2)/Mi 



dH 



= f{SK,SH), 



(A-7) 



where the function / only depends on the SK, and SH, 
because Mi and M2 only depends on the SIC and SH . By 
solving this differential equation, we obtain the magnetic 
field dependence of physical properties at zero tempera- 
ture. 
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